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Sperm swimming at low Reynolds number have strong hydrodynamic interactions when their concentration 
is high in vivo or near substrates in vitro. The beating tails not only propel the sperm through a fluid, but also 
create flow fields through which sperm interact with each other. We study the hydrodynamic interaction and 
cooperation of sperm embedded in a two-dimensional fluid by using a particle-based mesoscopic simulation 
method, multi-particle collision dynamics (MPC). We analyze the sperm behavior by investigating the relation- 
ship between the beating-phase diff'erence and the relative sperm position, as well as the energy consumption. 
Two eff'ects of hydrodynamic interaction are found, synchronization and attraction. With these hydrodynamic 
eff'ects, a multi-sperm system shows swarm behavior with a power-law dependence of the average cluster size 
on the width of the distribution of beating frequencies. 

PACS numbers: 82.70.-y, 87.16.Qp, 87.17.Aa 



I. INTRODUCTION 

Sperm motility is important for the reproduction of animals. 
A healthy mature sperm of a higher animal species usually has 
a flagellar tail, which beats in a roughly sinusoidal pattern and 
generates forces that drive fluid motion. At the same time, 
the dynamic shape of the elastic flagellum is influenced by 
the fluid dynamics. The snake-like motion of the tail propels 
the sperm through a fluid medium very efl&ciently. In the past 
decades, the effort to quantitatively describe the fluid dynam- 
ics of sperm has been very successful |[lll2]. 

However, despite considerable progress in modeling sperm 
elementary structures and the behavior of a single sperm in 
a fluid medium (H, 0], relatively few studies have examined 
the fluid-dynamics coupling of sperm and other mesoscopic or 
macroscopic objects, e.g., the synchrony of beating tails 
the tendency of accumulation near substrates H HI 111 
etc. In nature, the local density of sperm is sometimes ex- 
tremely high. For example, in mammalian reproduction, the 
average number of sperm per ejaculate is tens to hundreds 
of millions, so that the average distance between sperm is 
on the scale of ten micrometers — comparable to the length 
of their flagellum. The sperm are so close that the interac- 
tion between them is not negligible. In recent years, exper- 
iments JTll [13, EH H im O have revealed an interesting 
swarm behavior of sperm at high concentration, e.g. the dis- 
tinctive aggregations or 'trains' of hundreds of wood-mouse 
sperm 11141 [iSll . or the vortex arrays of swimming sea urchin 
sperms on a substrate |[T6ll . The mechanisms behind the abun- 
dant experimental phenomena are still unclear. In this paper, 
we focus on the hydrodynamic interaction between sperm and 
explain its importance for the cooperative behavior. 

The higher animal sperm typically have tails with a length 
of several tens of micrometers. At this length scale, viscous 
forces dominate over inertial forces. Thus, the swimming mo- 
tion of a sperm corresponds to the regime of low Reynolds 
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number [17]. Experimental observations of two Paramecium 
cells swimming at low Reynolds number have shown that 
the changes in direction of motion between two cells are in- 
duced mainly by hydrodynamic forces 1 18]. Studies of model 
micro-machines indicate that hydrodynamic interaction is sig- 
nificant when the separation distance is comparable to their 
typical size [19]. The hydrodynamic interaction between two 
rotating helices, like bacterial flagella, has been investigated 
both experimentally |(20l] and theoretically llillll]. An arti- 
ficial micro swimmer, which mimics the motion of a beating 
sperm, has been constructed from a red blood cell as head and 
a flagellum-like tail composed of chemically linked param- 
agnetic beads; the propulsion is then induced by a magnet- 
ically driven undulation of the tail ||23|] . Simulations have 
been employed to study the motion of a single of these arti- 
ficial microswimmer [24], as well as the hydrodynamic inter- 
actions between two swimmers Even studies of a mini- 
mal swimming model of three linearly connected spheres ||26|] 
have shown a complicated cooperative behavior i27|] . Thus, 
although there has been much progress on modeling and ob- 
serving a single swimmer, the understanding on the hydrody- 
namic coupling behavior of dense system of swimmers is still 
poor. 

In this paper, we focus on the cooperation behavior of 
sperm in two dimensions. Although real swimming sperma- 
tozoa are certainly three-dimensional, the qualitatively similar 
phenomena, and the great saving of simulation time, makes it 
worthwhile to discuss the problem of cooperation in a viscous 
fluid in two dimensions. Furthermore, sperm are attracted to 
substrates in in- vitro experiments iH HI, & 113] and are there- 
fore often swimming under quasi-two-dimensional condition 
(it has to be emphasized that hydrodynamic interactions in two 
dimensions and in three dimensions near a substrate are of 
course diflTerent). Thus, we construct a coarse-grained sperm 
model in two dimensions and describe the motion of the sur- 
rounding fluid by using a particle-based mesoscopic simula- 
tion method called multi-particle collision dynamics (MPC) 
||28[ 29] . This simulation method has been shown to capture 
the hydrodynamics and flow behavior of complex fluids over 
a wide range of Reynolds numbers very well IIO, EH], and is 



thus very suitable for the simulation of swimming sperm. 

This paper is organized as follows. Section II gives a brief 
description of our sperm model and of the particle-based hy- 
drodynamics approach. In order to understand a complex 
many-body system of micro- swimmers, a first-but-important 
step is to investigate the interaction between two swimming 
sperm. Thus, in Sec. Ill, we look at the cooperative behavior 
of two sperm. Two remarkable hydrodynamic eff'ects, syn- 
chronization and attraction, are found and discussed in detail. 
In Sec. IV, we analyze the clustering behavior of multi-sperm 
systems. In particular, we consider a sperm system with a 
distribution of beating frequencies, and determine the depen- 
dence of the cluster size on the variance of the frequency dis- 
tribution. 



II. SPERM MODEL AND MESOSCALE 
HYDRODYNAMICS 

A. Multi-Particle Collision Dynamics (MPC) 

MFC is a particle-based mesoscopic simulation technique 
to describe the complex fluid behaviors for a wide range of 
Reynolds numbers [SO, • The fluid is modeled by N point 
particles, which are characterized by their mass mi, continu- 
ous space position r/ and continuous velocity U/, where / = 
1,...N. In MPC simulations, time t is discrete. During every 
time step At, there are two simulation steps, streaming and 
collision. In the streaming step, the particles do not interact 
with each other, and move ballistically according to their ve- 
locities. 
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In the collision step, the particles are sorted into collision 
boxes of side length a according to their position, and inter- 
act with all other particles in same box through a multi-body 
collision. The collision step is defined by a rotation of all par- 
ticle velocities in a box in a co-moving frame with its center 
of mass. Thus, the velocity of the i-th particle in the j-th box 
after collision is 



Ui(t -h AO = Ucmjit) + - Ucmj] 



where 
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(3) 



is the center-of-mass velocity of j-th box, and is a ro- 

tation matrix which rotates a vector by an angle ±a, with the 
sign chosen randomly. This implies that during the collision 
each particle changes the magnitude and direction of its veloc- 
ity, but the total momentum and kinetic energy are conserved 
within every collision box. In order to ensure Galilean invari- 
ance, a random shift of the collision grid has to be performed 

The total kinematic viscosity v is the sum of two contri- 
butions, the kinetic viscosity Vkin and the collision viscosity 




FIG. 1 : (Color online) Two-dimensional model of sperm. The model 
consists of three parts, the head (blue), the mid-piece (red) and the 
tail (cyan). Two sinusoidal waves are present on the beating tail. 



Vcoii' In two dimension, approximate analytical expressions 
areillll. 
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where p is the average particle number in each box, m is the 
mass of solvent particle and h = At ^jhsT Ima^ is the rescaled 
mean free path. In this paper, we use ksT = l,m = I, a = 1, 
At = 0.025, a = 7i/2, p = 10. This implies, in particular, 
that the simulation time unit {ma^ jkBT^^ equals unity. With 
these parameters, the total kinematic viscosity of fluid is v = 
"^coii + "^kin ~ 3.02. The size of the simulation box is x Ly, 
where Lx = Ly = 200a, four times the length of the sperm tail, 
if not indicated otherwise. Periodic boundary conditions are 
employed. 



B. Sperm Model in Two Dimensions 

Although animal sperm difl'er from species to species, their 
basic structure is quite universal. Usually, a sperm consists 
of three parts: a head containing the genetic information, a 
beating long tail, and a mid-piece to connect head and tail. 
Our two-dimensional sperm model, shown in Fig.[Tl consists 
of these three parts. The head is constructed of Nhead = 25 
particles, where neighboring particles are linked by springs of 
finite length /q = 0.5a with interaction potential 
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into a circle of radius 2a. Each of the head particles has a 
mass rrihead = 20. The mid-piece consists of Nmid = 14 
particles of mass nimid = 10 connected by springs of length 
/o = 0.5a. The first particle of the mid-piece, which is fixed 
to the center of the head, is connected with every particle on 
the head by a spring of length I head-mid = 2a, in order to 
maintain the circular shape of the head, as well as to stabi- 
lize the connection between head and mid part. The tail has 
Ntaii = 100 particles of mass mtaii = 10, linked together by 
springs of length Iq. The spring constants are chosen to be 

head-mid = ^0"^ , head = , kmid = kail = 2 X 10^ whcrC 

head-mid IS the Spring constant for the connection of the head 
particles and the center, and k^id and head are the spring con- 
stants for the tail and the mid-piece, respectively. 

A bending elasticity is necessary for the mid-piece and the 
tail to maintain a smooth shape in a fluctuating environment. 
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and to implement the beating pattern. The bending energy is 
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where k denotes the bending rigidity, 9l(loCs) is a rotation ma- 
trix which rotates a vector anticlockwise by an angle foc^, and 
Cs,taii is the local spontaneous curvature of the tail of the ^- 
th sperm. We choose k = 10"^, much larger than the thermal 
energy ksT = 1 to guarantee that the mechanical forces dom- 
inate the thermal forces. For the mid part, the spontaneous 
curvature vanishes, Cs^mid = 0- <^s,taii ^ variable changing 
with time t and the position x along the flagellum to create a 
propagating bending wave. 



Cs,tail(^^ = Co,tail + ^ sin 
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A detailed analysis of the beating pattern of bull sperm in 
Ref . [361 shows that a single sine mode represents the beat to 
a very good approximation. The wave number q = Anll^Ntaii 
is chosen to mimic the tail shape of sea-urchin sperm [1], so 
that the phase difference between the first and the last parti- 
cles of the tail is An, and two waves are present (see Fig.[T]). fs 
is the beating frequency of the ^-th sperm. The constant c^^taii 
determines the average spontaneous curvature of the tail, ips is 
the initial phase of the first tail particle on the ^-th sperm, and 
A is a constant related to the beating amplitude. We choose 
A = 0.2, which induces a beating amplitude Ataii = 3.2a of 
the tail. As t increases, a wave propagates along the tail, push- 
ing the fluid backward at the same time propelling the sperm 
forward. We keep A, k, Ts, and (fs constant for each sperm 
during a simulation. Although the spontaneous local curva- 
ture is prescribed, the tail is elastic and its configuration is 
aff'ected by the viscous medium and the flow field generated 
by the motion of neighboring sperm. 

In order to avoid intersections or overlaps of diff'erent 
sperm, we employ a shifted, truncated Leonard- Jones poten- 
tial 
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between particles belonging to difl'erent sperm, where r is the 
distance between two particles. Parameters cr = I and e = 
13.75 are chosen. 

During the MPC streaming step, the equation of motion 
of the sperm particles is integrated by a velocity-Verlet algo- 
rithm, with a molecular-dynamics time step A^^ = 5 x 10"^, 
which is 1 / 50 of the MPC time step A^. The sperm only inter- 
acts with the fluid during the MPC collision step. This is done 
by sorting the sperm particles together with the fluid particles 
into the collision cells and rotating their velocities relative to 
the center-of-mass velocity of each cell. 

Since energy is injected into the system by the actively beat- 
ing tails, we employ a thermostat to keep the fluid temperature 



constant by reseating all fluid-particle velocities in a collision 
box relative to its center-of-mass velocity after each collision 
step. This procedure has the advantage that the energy con- 
sumption per unit time of the sperm can be easily extracted 
through the reseating of the particle velocities. 

We start with a single-sperm system with co,to// = and 
/ = 1/120. With the other parameters given in the previous 
section, our sperm model swims smoothly forward with the 
velocity Usingie = 0.016 ± 0.001. Because of its large size, 
the difl'usion coefficient of a sperm due to the thermal fluctu- 
ations of the MPC fluid is very small, on the order of 10"^ 
1I37I] . This implies that the time the sperm needs to cover a 
distance of half the length of its flagellum by passive difl'u- 
sion is more than a factor 10^ larger than the time to travel 
the same distance by active swimming. Therefore, diffusion 
plays a negligible role in our simulations. The energy con- 
sumption per unit time P single = 25.2 ± 2.4 — corresponding 
to P single - 3000kBTf. Thus, we estimate a Reynolds num- 
ber Re = 2Ataiiu single /y - 0-03 for our sperm model, where 
Ataii = 3.2a is the beating amplitude of the tail. 



III. TWO-SPERM SIMULATIONS 

A. Symmetric Sperm 

Two sperm, SI and S2, are placed inside the ffuid, initially 
with straight and parallel tails at a distance d = 5 (i.e. with 
touching heads). They start to beat at ^ = with different 
phases (fi and (f2. The initial positions of sperm do not matter 
too much, because two freely swimming sperm always have 
the chance to come close to each other after a sufficiently long 
simulation time. We consider two sperm with the same beat 
frequency / = 1/120, and the same spontaneous curvature 

C0,tail - ^^ 

In the dynamical behavior of these hydrodynamically inter- 
acting sperm, two eflfects can be distinguished, a short time 
"synchronization" and a longer time "attraction" process. If 
the initial phase difference A^? = ^2 - at time ^ = is 
not too large, an interesting effect denoted "synchronization" 
takes place, which is accomplished within a few beats. This 
process is illustrated in Fig.|2t-e by snapshots at different sim- 
ulation times. The synchronization time depends on the phase 
difference, and varies from about two beats for A^? = 0.57r 
(see Fig. [2]) to about five beats for A(f = n. A difference 
in swimming velocities adjusts the relative positions of the 
sperm. After a rapid transition, the velocities of two cells be- 
come identical once their flagella beat in phase. Because the 
initial distance between tails d = 5 is smaller than the beating 
amplitude 2Ataii = 6.4a, the sperm tails can touch when they 
start to beat for O.671 < Acp < I. An. This geometrical effect 
is reduced by the hydrodynamic interaction, which aff'ects the 
beating amplitude. In case contact occurs, it accelerates the 
synchronization. In order to avoid this direct interaction due to 
volume exclusion, we have also performed simulations of two 
sperm with initial distance J = 10, and find the synchronized 
state achieved within several beats, as in the simulations with 
d = 5. Thus, the synchronization effect is of purely hydrody- 



4 




(g) 



FIG. 2: (Color online) Snapshots of two sperm with phases (p\ (up- 
per), if2 (lower), and phase difference = ^2 - V^i - 0.5:7r. (a) 
tf = 1/6 initial position; (b) tf = 2/3; (c) tf = 7/6; (d) tf = 5/3; 
(e) tf = 13/6; (f) tf = 10^; (g) tf = 100^. From (a) to (e), the 
synchronization process takes place. The tails are already beating in 
phase in (e). From (e) to (g), two synchronized sperm form a tight 
cluster due to hydrodynamic attraction. 



namic origin. Since the beating phase at time t is determined 
by / and (ps, which are kept constant in our simulations, our 
model sperm can only achieve synchronization by adjusting 
the relative position. 

Our results are in good agreement with the prediction of 
Taylor f3], based on an analytical analysis of two-dimensional 
hydrodynamics, that the viscous stress between sinusoidally 
beating tails tends to force the two waves into phase. The 
same phenomenon has also been observed by Fauci and Mc- 
Donald 10] in their simulations of sperm in the presence of 
boundaries, and has been called "phase-locking" effect. A 
similar effect of undulating filaments immersed in a two- 



dimensional fluid at low Reynolds number was seen by Fauci 
in Ref. |38]. 

Synchronization is a fast process, which is achieved in at 
most ten beats in our simulations. Another hydrodynamic ef- 
fect, which we denote "attraction", takes much longer time. 
Two synchronized and separated sperm gradually approach 
each other when they are swimming together, as if there was 
some effective attractive interaction between them. The only 
way in which the sperm can attract each other in our simula- 
tions is through the hydrodynamics of the solvent. This ef- 
fect takes several ten beats to overcome the initial distance of 
d = 5 between the tails. The final state of attraction, in which 
the sperm tails are touching tightly, is shown in Fig. [2^. 

Fauci and McDonald did not see the hydrodynamic at- 
traction, because they considered a sperm pair confined be- 
tween two walls. As explained in Ref. |6|, there is an evident 
tendency for a single sperm to approach the wall. When two 
parallel sperm are placed between the walls, there seems to 
be a critical initial distance between the sperm, below which 
synchronization occurs, and above which swimming towards 
the wall occurs. Our understanding is that, in their simula- 
tions, the viscous drag towards the walls was competing with 
the viscous attractive effect between sperm. Hence in some 
cases, they could only see a synchronization effect, and nei- 
ther a clear towards-wall tendency nor a distinguishable at- 
traction effect. The hydrodynamic attraction was masked by 
the presence of the walls. 

To analyze the cooperating sperm pair in more detail, we 
choose the head-head distance dh to characterize the attrac- 
tion and synchronization, because it is easy experimentally to 
track the head position. The dependence of dh on the phase 
difference is symmetric with respect to A(fi = because of the 
symmetry of the sperm structure. Thus, we show in Fig. [3] 
only results for A(f > 0. There is a plateau at about dh = 5a 
for A(f < OAn, which corresponds to the sperm heads touch- 
ing each other. For A(p > 0.4;r, dh increases linearly with A(fi. 
Finally, for A(f > l.Sn, the phase difference is so large that 
the attraction is not strong enough to overcome the thermal 
fluctuations and pull the sperm close together. Although syn- 
chronization still occurs at the beginning, the two sperm leave 
each other soon after. 

Riedel et al. (l^ also see such a linear relation in their ex- 
periments of sea-urchin sperm vortices. They define the beat- 
ing phase of a sperm by its head oscillation, and an angular 
position of the sperm head within the vortex. In this way, the 
beating phase difference of the sperm in the same vortex was 
found to have a linear relation with the angular position dif- 
ference, which corresponds to the head-head distance in our 
simulations. 

So far, we have considered sperm with a single beat fre- 
quency. In nature, sperm of the same species always have a 
wide distribution of beat frequencies. For example, the beat 
frequency of sea-urchin sperm ranges from 30Hz to 80Hz 
ill], and the frequency of bull sperm ranges from 20Hz to 
30Hz ||36|] . Thus, we assign different beat frequencies to 
two sperm, fi = 1/120 and /2 = 1/119.4, corresponding to 
A///1 ~ 0.5%, but set the same initial phases (fs = 0. This 
implies that the phase difference of the beats between the two 
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FIG. 3: (Color online) Head-head distance dh of two cooperating 
sperm. Simulation data are shown for fixed phase diff'erence (red, 
□), with error bars denoting the standard deviation. The interpolating 
(red) line is a linear fit for OAn < Aip < l.Sn. The distance dh is also 
shown as a function of time t (top axis) in a simulation with a 0.5% 
diff'erence in the beat frequencies of the two sperm (solid line). 



FIG. 4: Snapshots of two synchronized human sperm in experiment 
at diff'erent times |[3^ . |4Q|i . (Left) Two sperm with initially well syn- 
chronized tails and very small phase diff'erence; (Middle) the sperm 
are still swimming together and are well synchronized after 4 sec- 
onds; a phase diff'erence has developed; (Right) the sperm begin to 
depart after 7 seconds. The scale bar corresponds to a length of 25 



sperm increases linearly in time, 

Aif = 27l(f2-fl)t. (11) 

Fig. [3] shows the head-head distance versus time. It agrees 
very well with the data for fixed phase differences. At tfi = 
150 where A(fi ^ 1.5;:, the sperm trajectories begin to depart. 

Fig. [4] shows two cooperating human sperm swimming in 
an in- vitro experiment near a glass substrate lHH H^]. The 
two sperm swim together for more than 6 seconds at a beat 
frequency of approximately 8Hz. Their tails remain synchro- 
nized during this time, while the head-head distance and phase 
difference increases with time (see Fig.©. After a while, the 
sperm leave each other because the phase difference becomes 
too large. There is no indication of a direct adhesive interac- 
tion between the sperm. 

An interesting question is whether the cooperation of a 
sperm pair reduces the energy consumption. Fig. [5] displays 
the energy consumption of two sperm with the same beat fre- 
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FIG. 5: (Color online) Energy consumption per unit time, P, of two 
cooperating sperm. Symbols show simulation data for fixed phase 
diff'erence (red, □), where error bars denote the standard deviation. P 
versus time nn a simulation with a 0.5% diff'erence in the frequencies 
of two sperm (solid black line). 



quency /= 1/120 as a function of the phase difference. The 
energy consumption P is nearly constant at small phase dif- 
ference. It increases for A(p > O.Sn roughly linearly until it 
reaches another plateau for A(f > l.Sn. The second plateau 
corresponds to two sperm swim separately, so that energy con- 
sumption is twice the value of a single sperm. Our results are 
in agreement with the conclusion of Taylor ^ that less energy 
is dissipated in the fluid if the tails are synchronized. 

Fig. [5] also shows the energy consumption of two sperm 
with fi = 1/120, /2 = 1/119.4 and = ^2 = as a function 
of time t. In this simulation, we start with two sperm which are 
parallel and at a distance d = 5. For tfi < 25, the energy con- 
sumption decreases as the sperm are approaching each other. 
The data agrees quantitatively very well with results for con- 
stant A(fi, and reaches a plateau when the cooperating sperm 
pair departs. 

The synchronization and attraction also exists in our sim- 
ulation of swimming flagella without heads. In this case, the 
time-reversal symmetry of Stokes flow implies that no syn- 
chronization nor attraction is possible at zero Reynolds num- 
ber. In our simulations, the thermal fluctuations and a finite 
Reynolds number break the time-reversal symmetry. 



B. Asymmetric Sperm 

In nature, sperm have an abundance of different shapes. In 
particular, these shapes are typically not perfectly symmet- 
ric. The asymmetric shape can cause a curvature of the sperm 
trajectory [7, 41]. For example, sea-urchin sperm uses the 
spontaneous curvature of the tail to actively regulate the sperm 
trajectory for chemotaxis ||42[|43|] . In our simulations, we im- 
pose an asymmetry of the tail by employing a non-zero spon- 
taneous curvature co^taii- 

We consider curved sperm tails, with cojaii = 0.04/a, which 
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FIG. 6: (Color online) Head-head distance dh of two cooperating 
sperm with spontaneous curvature co,ra// - 0.04/a as a function of the 
phase difference A^. The error bars represent standard deviations. 
Lines are linear fits to the data in the range -\.97i < I^(p < -Q.5n 
and 0.3:7r < A(p < O.Itt, respectively. The inset shows two typical 
conformations with positive and negative phase diff'erence. 
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FIG. 7: (Color online) Energy consumption per unit time, P, versus 
phase diff'erence Aip of two sperm with spontaneous curvature co^taii - 
0.04/a of their tails. The error bars represent standard deviations. 



results in a mean curvature of the trajectory of a single sperm 
oi CfU = 0.041 ± 0.009. For sperm with curved tails, the 
head-head distance dh(A(f) is not symmetric about A(f = 0, 
as shown in Fig. [6l Here A(f is defined as the phase of the 
sperm on the inner circle minus the phase of the sperm on the 
outer circle. The steric repulsion of the heads causes a plateau 
of the head-head distance at dh = 5 for small phase dififer- 
ences A^, as in Fig. [3] for symmetric sperm. For A(fi < -nil 
and Aip > 7r/4, the head distance increases linearly with in- 
creasing phase difiference, with a substantial difiference of the 
slopes for A^ < and Aip > 0, see Fig. [6l The two sperm 
depart when Aip > O.ln. For Aip < -l.On, the sperm pair 
briefly looses synchronicity, but then rejoins with a new phase 



diflTerence Aip' = Aip -\- In. 

The energy consumption P for sperm with spontaneous cur- 
vature (see Fig. [7]) also increases sharply at Aip = O.Stt and 
stays at the plateau with P = 51.0 ± 2.8 for larger Aip, as the 
sperm are swimming separately. However, for Aip < -7i/2, P 
increases rather smoothly until the cooperation is lost for large 
phase diflTerences. 

We conclude that the strong curvature of the tail breaks the 
symmetry of the head-head distance dh and the energy con- 
sumption P in Aip, but the eflTect of synchronization and at- 
traction are still present and play an important role in the co- 
operation of sperm pairs. 

IV. MULTI-SPERM SYSTEMS 

When two sperm with the same beating period happen to 
get close and parallel, they interact strongly through hydro- 
dynamics and swim together. With this knowledge of hydro- 
dynamic interaction between two sperm, we now study a sys- 
tem of 50 sperm in a simulation box of 200 x 200 collision 
boxes, which corresponds to a density of about three sperm 
per squared sperm length. The initial position and orienta- 
tion for each sperm are chosen randomly. Considering that in 
real biological systems the beat frequency is not necessary the 
same for all sperm, we perform simulations with Gaussian- 
distributed beating frequencies. The initial phases of all sperm 
are ips = 0. 

We consider a system of symmetric sperm. Fig. [8] 
shows some snapshots of systems with diflTerent width Sf = 
{(Af)^y^^/{f) of the Gaussian frequency distribution. Here, 
((A/)^) is the mean square deviation of the frequency distri- 
bution, and (/) = 1/120 is the average frequency. 

For Sf = 0, once a cluster has formed, it does not disinte- 
grate without a strong external force. A possible way of break- 
up is by bumping head-on into another cluster. For Sf > 0, 
however, sperm cells can leave a cluster after sufi&ciently long 
time, since the phase diflTerence to other cells in the cluster in- 
creases in time due to the dififerent beat frequencies (compare 
Sec. lIII Al) . At the same time, the cluster size can grow by col- 
lecting nearby free sperm or by merging with other clusters. 
Thus, there is a balance between cluster formation and break- 
up, as shown in the accompanying movie ll44l] . Obviously, the 
average cluster size is smaller for large Sf than for small Sf 
(see Fig. [8]). 

To analyze the multi-sperm systems, we define a cluster as 
follows. If the angle between vectors from the last to the first 
bead of the tails of two sperm is smaller than 7r/6, and at the 
same time the nearest distance between the tails is smaller 
than 4a, which is approximately 1/10 of the length of the tail, 
then we consider these two sperm to be in the same cluster. 
By this definition, we find the evolution of the average cluster 
size < He > shown in Fig. |9l Here, < > is the average 
number of sperm in a cluster, 

< ric >=^ncU(nc) , (12) 

ric 

where Yl(nc) is the (normalized) cluster-size distribution. For 



7 




FIG. 8: (Color online) Snapshots from simulations of 50 symmetric 
sperm with different widths (5/ of a Gaussian distribution of beating 
frequencies, (a) 6f = 0; (b) 6f = 0.9%; (c) 6f = 4.5%. The red 
ellipses in (b) and (c) indicate large sperm clusters. The black frames 
show the simulation boxes. Note that we employ periodic boundary 
conditions. 

Sf = 0, the average cluster size continues to increase with 
time. Both systems in Fig. [9] with ^/ > reach a stationary 
cluster size after about 50 beats. The stationary cluster size is 
plotted in Fig.[TOlas a function of the width ^/ of the frequency 
distribution. We find a decay with a power law, 

< He >~ sy (13) 

with 7 = 0.20 ±0.01. The error for y is estimated from a fit of 
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FIG. 9: (Color online) Time dependence of the average cluster size, 
< ric >, in a system of 50 symmetric sperm with various widths ^/ of 
the frequency distribution, as indicated. 
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FIG. 10: (Color online) Dependence of the average stationary cluster 
size, < ric >, on the width of the frequency distribution 6f. Data are 
shown for a 50-sperm system (■) and a 25-sperm system (o). The 
lines indicate the power-law decays <nc>- 2.12(5"^^^^ (upper) and 
<nc >= 1.55(5-0-196 (lower). 



the data for both 50-sperm and 25-sperm systems. The neg- 
ative power law indicates that the cluster size diverges when 
^/ ^ 0. This tendency is also implied by the continuously 
increasing cluster size for = in Fig.|9l The cluster- size 
distribution in the stationary state is shown in Fig.[TT] Cluster- 
size distributions have been studied in much simpler systems 
of self-propelled particles, such as point particles with a con- 
stant magnitude of their velocities, which adjust their travel- 
ing direction to the direction of their neighbors | 45]. In such 
simplified models, a power-law decay of the cluster- size dis- 
tribution with an exponent in the range of 1.5 to 1.9 have 
been found |46], followed by a rapid decay for large clus- 
ter sizes due to finite-size eff'ects. Similarly, in a system of 
self-propelled rods with volume exclusion, a crossover from 
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FIG. 11: (Color online) Cluster size distribution, Il{nc). Data are 
shown for 25 sperm in a 200 x 200a^ box (•), 100 sperm in a 400 x 
400a^ box (□) [note that both systems have the same sperm density], 
and 50 sperm in a 200 X 200a^ box (A). The lines correspond to an 
exponential distribution. The inset shows the same data in a double- 
logarithmic representation. The line indicates a power law n~^-^. 

power-law behavior at small cluster- sizes to a more rapid de- 
cay for large cluster sizes has also been found 1470 . A power- 
law decay of the cluster- size distribution is indeed consistent 
with our results for smaller cluster sizes, as shown in the inset 
of Fig. [TT] The rather similar value of the exponent with that 
of Ref. ll46ll is probably fortuitous. We attribute the exponen- 
tial decay of the cluster-size distribution for larger cluster size, 
which is apparent in Fig.[TTJ to finite-size eff'ects. Simulations 
of larger system sizes are required to confirm this conclusion. 

To analyze the energy consumption of sperm clusters, we 
consider a special case where sperm of the same frequency are 
prearranged to pack tightly and to be synchronized, as shown 
in Fig. [T2h . A simple linear relationship between the energy 
consumption of the sperm cluster and the cluster size is shown 
in Fig.[T2h. From the linear fit of the data, we obtain an energy 
consumption per sperm for an infinitely large cluster, P/ric = 
13.7. Thus, a freely swimming sperm can reduce its energy 
consumption by almost a factor 2 by joining a cluster. 

The swimming speed of a sperm cluster decreases slowly 
with increasing sperm number, as shown in Fig. [T2b . When 
flagella are very close, with distances smaller than the size of a 
MPC collision box, hydrodynamic interactions are no longer 
properly resolved. Instead, the collision procedure yields a 
sliding friction for the relative motion of neighboring flagella. 
Thus, the energy of the beat is not only used for propulsion, 
but also to overcome the sliding friction. The energy con- 
sumption of tail-tail friction is proportional to the number of 
neighbor pairs, and the hydrodynamic resistance of moving 
the whole cluster is proportional to the cluster size and speed. 
Thus, the total energy consumption can be written as 

P = Cny^Pf(nc-l\ (14) 

where Pf is the energy consumption due to tail-tail friction, 
and C is a constant. With the relation P = 13.7/2^ + 10.1 



30O 




c 



FIG. 12: (Color online) (a) Energy consumption per unit time, P, 
of sperm clusters as a function of cluster size ric. Symbols indicate 
simulations results. The fit line (red) is given by P = 13.7/Zc + 10.1. 
The inset shows an illustration of an arranged cluster of 20 sperm, 
(b) Center of mass speed of sperm clusters as a function of cluster 
size He. The fit line (red) is given by v = 0.00334 V6.42 + 16.4/ ric. 



obtained above, the data for the cluster speed can be fitted 
to Eq. (El), which yields pf = 12^ and C = 8.96 x lOl 
Thus, the cluster speed reaches a non-zero asymptotic value 
[(13.7 - Pf)/CY^^ ^ 0.0085. for large cluster size, about a 
factor 2 slower than a single sperm. 



V. SUMMARY AND DISCUSSION 

We have simulated the hydrodynamic interaction between 
sperm in two dimensions by the multi-particle collision dy- 
namics (MPC) method. Two eff'ects of the hydrodynamic 
interaction were found in our simulations. First, when two 
sperm are close in space and swimming parallel, they syn- 
chronize their tail beats by adjusting their relative position. 
This process can be accomplished in a very short time, less 
than 10 beats. Second, two synchronized sperm have a ten- 
dency to get close and form a tight pair. This process takes 
much longer time then synchronization. It usually takes about 
100 beats to overcome a distance of 1/10 tail length between 
sperm in our simulations. 



9 



These hydrodynamic effects favor the cooperation of sperm 
in motile clusters. For a multi-sperm system, the average clus- 
ter size diverges if all sperm have the same beating frequency. 
A distribution of frequencies leads to a stationary cluster- size 
distribution with a finite average cluster size, which decreases 
with a power law of the variance of the frequency distribution. 
Furthermore, the average cluster size increases with increas- 
ing sperm density. The probability to find a cluster decreases 
with a power law for small cluster sizes; an exponential decay 
for large cluster sizes is attributed to finite- size eff'ects. 

In sperm experiments, large bundles have been found in 
some species, like fish flies IllKlsl] and wood mouse |[l2[[l5|] . 
For fish-fly sperm, this has been attributed to some agglu- 
tination of the sperm heads to keep the size and structure 
of the bundles. Wood mouse sperm were released into an 
in-vitro laboratory medium, initially in single cell suspen- 
sion 1IT2I] . Within 10 minutes, large bundles containing hun- 
dreds or thousands of sperm were formed as motile 'trains' 
of sperm. Motile bundles of 50-200 sperm were also found 
in the after-mating female's body, as well as many non-motile 
single sperm. The hook structure on the head of wood mouse 



sperm is believed to favor the formation of such huge cluster 
in in-vitro experiments. 

In our simulations, sperm clusters are always seen, e.g. 
as marked in Fig. [H after the system has reached a dy- 
namically balanced state of cluster sizes. Thus, we predict 
that hydrodynamic synchronization and attraction play an 
important role in the cluster formation of healthy and motile 
sperm, such as the bundles and trains observed for fish-fly 
and wood-mouse sperm at high concentrations, respectively. 
Furthermore, since the cluster size decreases with increasing 
width ^/ of the distribution of beat frequencies, our results 
are consistent with the experimental observation that if the 
sperm are hyperactivated 1IT2I1 . which is an abnormal beat 
mode, or if some sperm are dead, the clusters fall apart. 
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